Nonuniqueness in spin-density-functional theory on lattices 
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In electronic many-particle systems, the mapping between densities and spin magnetizations, 
{n(r),m(r)}, and potentials and magnetic fields, {v(r), B(r)}, is known to be nonunique, which 
has fundamental and practical implications for spin-density- functional theory (SDFT). This paper 
studies the nonuniqueness (NU) in SDFT on arbitrary lattices. Two new, non-trivial cases are 
discovered, here called local saturation and global noncollinear NU , and their properties are discussed 
and illustrated. In the continuum limit, only some well-known special cases of NU survive. 
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Density-functional theory (DFT) 0, is a widely 
used approach for calculating the electronic structure of 
atoms, molecules, and many types of materials. The fun- 
damental theorem of Hohenberg and Kohn 0] establishes 
that the ground-state density n(r) of a nonmagnetic, 
nondegenerate electronic system uniquely determines the 
scalar potential v(r), apart from an arbitrary additive 
constant C. However, many applications of interest in- 
volve electronic ground states that have a spin magne- 
tization, m(r), and/or are exposed to external magnetic 
fields, B(r). Such situations can be handled with spin- 
density-functional theory (SDFT) 0, H 0, @ , where 
the fundamental variable is the 4-density {n(r), m(r)}, 
which couples to the 4-potential {v(r), B(r)}. 

Using the Rayleigh-Ritz variational principle, one can 
prove that a given physical 4-density uniquely determines 
the ground-state wave function Wo (apart from trivial 
phase factors). This guarantees that any property of the 
system expressible in terms of Wo is a functional of the 4- 
density. Following the Hohenberg-Kohn theorem of DFT, 
one might expect that there is also a unique map from 
ground-state wave functions to 4-potentials. However, 
it was recognized long ago by von Barth and Hedin Q, 
and more recently by Capelle and Vignale [f| and Esch- 
rig and Pickett [y] that such a unique correspondence 
does not exist. These authors showed that the extent of 
nonuniqueness (NU) in SDFT is much greater than the 
addition of a mere constant C to v(r) for the nonmag- 
netic systems considered in DFT. As a consequence, some 
applications of SDFT such as the calculation of excita- 
tion energies or of one-electron spin gaps in half-metallic 
ferromagnets must be critically reexamined. 

NU in SDFT means that an TV-particle ground state 
Wo remains unchanged under addition of a 4-potential 
(Aw(r), AB(r)}. This happens if and only if Wo is an 
eigenstate of the operator yj 



JV 



AH = Y,[Av(r j )-AB(v j ).a j ] 



(1) 



where <jj is the vector of Pauli matrices acting on the spin 
of the jth electron, and we use units where the Bohr mag- 



neton fiB = 1. The entire spectrum {W} is invariant if 
AH is related to a constant of motion (systematic NU); 
all other cases are called accidental NU |fj . Similar con- 
ditions for NU can be given for other multi-component 
generalizations of DFT, e.g. for current-carrying or su- 
perconducting systems |l0|. Thus, NU appears to be a 
characteristic feature of generalized DFTs, and its fun- 
damental and practical consequences need to be explored 
in detail. In the first place, it is important to know what 
types of NU can occur in practice. The following classes 
of examples have been identified in SDFT: 

(a) {Av = C, AB = Be z }, with constant B, for sys- 
tems with an energy gap and a collinear spin arrangement 
along the z-axis, where the {W} are eigenstates of S z (sys- 
tematic NU). B should be sufficiently small to avoid level 
crossings in order for Wnto remain the ground state. As 
shown by Gidopoulos the mapping in the collinear 
case is unique in the broader sense that spin-potentials 
{u-f (r), V[ (r)} which differ by more than a spin-dependent 
constant always have different ground states. 

(b) In a fully spin polarized Kohn-Sham system with 
n|(r) = n(r) and nj_ = and an energy gap, there is an 
infinite number of spin-down Kohn-Sham potentials that 
produce the same ground state (accidental NU). 

(c) {Av — Au(r),AB = Au(r)m(r)/m(r)} for one- 
electron systems only, where it(r) is an arbitrary function 
and A is sufficiently small (accidental NU) |4|. 

In Ref. [||, a general condition for NU in A'-electron 
systems is given: If Wo is invariant under the addition of 
a 4-potential, then this 4-potential must have the form 
{Av(r) = C, AB(r) = Be(r)}, i.e., AB may possibly be 
noncollinear but must be constant in magnitude. How- 
ever, Argaman and Makov [7( have raised doubts whether 
such noncollinear 4-potentials can really be found. 

It thus appears that NU for N > 1 electrons, while in- 
teresting and of potential practical relevance, is limited 
to rather simple situations of collinear spin arrangement 
or full, ferromagnetic spin polarization. The purpose of 
this paper is to examine NU in SDFT for arbitrary lat- 
tice systems. This has technical advantages over dealing 
with continuum systems since one can use linear algebra 
methods in finite vector spaces. We discover two new, 
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non-trivial classes of NU for iV-electron systems with 
noncollinear spins, both of the "accidental" type. These 
examples require the ground state 4-density to satisfy cer- 
tain constraints on the lattice. In the continuum limit, 
we show that only some well-known special cases survive, 
and we discuss consequences for practical applications. 

We consider a noninteracting TV-electron system on 
a finite-size lattice with P lattice points whose specific 
geometry is not important for the following. We as- 
sume that the kinetic-energy operator T has been suit- 
ably discretized on this lattice, for example using a finite- 
difference approach. The single-particle wave functions 
tpj obey the following Schrodinger equation: 



[f + V-B- a]tpj = Ejtpj 



J 



1. 



,2P. 



(2) 



For the spatial part of the wave functions we use a local- 
ized basis, (pug = Siki i, k = 1, . . . , P. The jth eigenstate 
on lattice site k can then be written as 



jk 



P 
i=l 



pot + djiP) 



c jk a + d jk P 



(3) 



where a, /3 are the usual two-component spinors. The 
coefficients Cj k ,dj k follow from diagonalizing the Hamil- 
tonian matrix associated with the lattice-specific T and 
the given 4-potential {w,B}. The resulting 4-density on 
lattice site k is 



( n k \j 

m 
m 



W/ 



N 

£ 

3=1 



/ \c ]k \ 2 + \d ]k \ 2 ^ 
Cjkdjk + c* jk d jk 



ic jk d* k 



ic- k d jk 



(4) 



\ \c jk \ 2 - \d jk \ 2 J 



A trivial case of NU arises when the lattice holds the 
maximum number of electrons allowed by the Pauli prin- 
ciple, N = 2P, for then n = 2 and m = for all external 
4-potentials. We now formulate the first new nontrivial 
example for NU on lattice systems, which we call local 
saturation NU. 

Theorem I: A noninteracting N -electron ground state 
on a P -point lattice is invariant under a perturbation with 
4-potential {v' kl B' k = v k m k /m k }, with arbitrary v' k , that 
acts locally only on those lattice sites where n k = m k . 

Notice that for N = 1 this reduces to the example (c) 
above, since in that case n = m on all lattice sites. For 



N > 1, it is not hard to show that the condition n k = m k 
on a specific site k is satisfied if and only if 



1>jk = fj(C k a + D k 0), j = l,. 



(5) 



i.e., the lowest N ipj must have the same spin part, and 
therefore parallel magnetizations, on point k. Let us now 
act on these states with a 4-potential that is nonzero on 
site k only and vanishes on all other sites i =/= k: 

{v' k -B' k -S)f 3 {C k a + D k (3) 

- I 3 [C k v' k - D k B' x ^ k + iD k B' y<k - C k B'^ k ]a 

+ j:ihr[ ( */. /''.. /, - iC k B' yM + l),ir : X i . (6) 
From this, it is straightforward to show that 

v' k [l - (m k /n k ) ■ a]f 3 {C k a + D k f3) = (7) 

for each j, which proves Theorem I. Local saturation NU 
can occur on one or more isolated lattice sites, but also on 
groups of sites, which includes examples associated with 
the formation of ferromagnetic domains on the lattice. 

The second new class of examples belongs to the 
Eschrig-Pickett type [§j and will be referred to in the 
following as global noncollinear NU. 

Theorem II: A noninteracting 2- electron ground state 
on a P -point lattice is invariant under a perturbation with 
4-potential {v' = 0,B' = Am/m}. |A| = const., if the 
ground state satisfies sign(\)(fii — n-i)lm — const. 

Here, fi\ and n% denote the two occupied orbital den- 



sities, with n 



ik 



\Cjk\ 



1,2. To prove 



Theorem II, one needs to show that the magnetic field 
B' = Am/m causes at most an orthogonal rotation 
within the space spanned by the two lowest single-particle 
eigenstates, ipi and -02, which leaves the associated 2- 
particle Slater determinant invariant. Thus, 



H a -X- 



(8) 

for % = 1, 2, where the 7™ form an orthogonal 2x2 ma- 
trix, and Ho is the unperturbed single-particle Hamilto- 
nian whose first two eigenstates and energies are 12 and 
Ei 2, see Eq. (J5J. Now consider a lattice site k. After 
some straightforward algebra, using rri^ — m x ±im v , one 
arrives at the following expression: 



cikEi 


_ ht_ 






m k 






d lk Ex 






m k 



(cifcWfc + d lk m k ) 



C2fc^2 - — {c 2k m z k + d 2k m k 
m k 

d 2k E 2 -(c 2k m~l - d 2k m z k 

m k 




(9) 



which has the form of a generalized 2x2 eigenvalue problem. Eq. © can be easily transformed into a standard 
eigenvalue problem by multiplying with the inverse of the right-hand coefficient matrix. Using relations |@J , one finds 
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after some manipulation 



(nik 

m k 

A A 



ri2k) 



- 2 {ci k c* 2k + d lk d* 2k ) 



2 ( c ife c 2fc + d* lk d 2k ) 

m k 

— {n\k ~ n 2 k) + E 2 
m k 







n 




\ l2i ) 





(10) 



which leads to a characteristic second-degree polynomial 
with solution 

Ei+E 2 i{AEf .„n lk -n 2k 
= — 5— ± ] l'—+Xl-X k AE— — , 

(11) 

where AE = E\—E 2 . We see immediately that for A = 
this reduces to £1.2 = -Ei,2- 

So far, the derivation was for a specific lattice site k. 
To ensure that the solution e 12 and the associated or- 
thogonal eigenvectors 7^ are the same for all P lattice 
sites, we need to impose the constraints 

. .fhn-fhi . ,nip-n 2P 

sign(Ai) = . . . = sign(Ap) , (12) 

mi mp 

which completes the proof of Theorem II, and deter- 
mines sign(Afc). Global noncollinear NU thus requires 
two-electron ground states whose orbital densities and 
total magnetization are related according to Eq. (|12[) . 
An explicit example for this will be given below. 

Furthermore, from the normalization of the orbital 
densities, 2i=i(^ii ~^2i) = 0, one finds that global non- 
collinear NU requires a total magnetization of the form 

sign(Ai)mi = -sign(A 2 )m 2 - sign(A 3 )m 3 

— ... — sign(Ap)mp . (13) 

Some additional remarks are in order: 

One-electron case. For a single electron, |j5Jl reduces to 

e = Ei- \ k n k /m k . (14) 

Again, we require this to be the same on all lattice sites. 
But, of course, n = m everywhere for a single electron, 
so that we end up with the condition A = const., i.e. 
e = Ei — A. This leads to the statement that any one- 
electron ground state is unchanged under the influence 
of a magnetic field B' = Am/m (provided A is suffi- 
ciently small such that the order of the lowest levels is 
not changed) . This is a special case already contained in 
the one-electron limit of Theorem I. 

N -electron case. Global noncollinear NU cannot occur 
for systems with more than two electrons, which can be 
seen as follows. The 2-electron derivation is easily gener- 
alized up to the point where one arrives at a generalized 
eigenvalue problem similar to Eq. ©, but of the type 
Rrj = eSfj where R,S_ are 2 x N rectangular matrices, 



and 7 is an ./V-component column vector. Such under- 
determined problems are singular, that is, one can find 
at most two ei gen values, all remaining N — 2 eigenvalues 
are undefined [llj. Often one finds no eigenvalues at all. 
This means that, except for trivial situations or by acci- 
dent, there is no noncollinear field B' that results only in 
a rotation within the single-particle ground-eigenspace. 
The iV-particle ground-state Slater determinant is thus 
not invariant for N > 2. 

2-point lattices. For the special case P = 2 one can 
show that all well-behaved 4-potentials produce two- 
electron ground states whose magnetization has same 
magnitude on the two lattice sites, m\ — m 2 . This re- 
sult is independent of electron interactions. Global non- 
collinear NU is thus always present on 2-point lattices. 

An interesting implication of this is that magnetiza- 
tions with mi ^ m 2 can arise on a 2-point lattice only as 
ensemble 4-densities of degenerate ground states. Similar 
consequences of NU are expected for lattices with more 
than 2 points, i.e., certain classes of 4-densities can only 
come from ensembles of degenerate ground states. For 
the case of non-magnetic DFT, the topology of the v- 
and n-spaces on lattices was recently clarified 13], with 
the result that pure- and ensemble-u-representable den- 
sities have the same mathematical measure. In SDFT, 
this general statement no longer holds due to the much 
richer NU, as is evident from the 2-point lattice example. 

Fig. ^ illustrates an example for global noncollinear 
NU on a linear 3-point lattice with lattice constant a 
and sites 1,2,3, using a finite-difference T. To discover 
this and many other examples, the 4-potential parame- 
ter space was numericall y s earched with a multidimen- 
sional simplex algorithm [l2J until a two-electron ground 
state was found to satisfy Eq. 112f) . with mi = m 2 + 7713, 
to within an accuracy of 1CU 14 (similar numerical tech- 
niques yield examples for local saturation NU). Measur- 
ing energies in units of Ti 2 /2mo? , we give the 4-potential 
and resulting 4-density in Table HI All magnetic fields 
B' = B ± Am/rn (+ on site 1, — on sites 2,3) produce 
the same two-electron ground state 4-density (keeping v 
fixed), for —1.7 < A < 0.7. Values of A outside that range 
result in level crossings and thus different ground states. 

We now turn to the continuum limits of our lattice 
examples. 

1. Local saturation NU. It is possible that the local 
condition n(r) = m(r) is satisfied in lower-dimensional 
subspaces (e.g., points or lines) for a continuum system. 
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FIG. 1: Global noncollinear NU for a linear 3-point lattice 
with mi = rri2 + m3. All magnetic fields shown here, where 
= Bi+Ami/mj and B 2 . 3 = B2,3 — Am2,3/m2,3 (see Table 
Pi, produce the same 4-density (keeping v fixed). 

TABLE I: The 4-potential and 4-density used in Fig. Q] 
on lattice sites 1, 2, and 3. Potentials and magnetic fields 
are measured in units of h 2 /2ma 2 (setting /ib = 1). The 
lattice 4-density is dimensionless. Constraints If 21 and If 311 
are satisfied, with |ni — n^jm — 0.99756 on each lattice site. 





1 


2 


3 


V 


-1.62192 


1.55381 


0.0 


B x 


0.87156 


-0.14000 


-0.50808 


B y 


0.15523 


0.23990 


-0.69702 


B z 


1.76179 


-0.77994 


0.11107 


ni,n 2 


0.96869,0.07157 


0.02926,0.20498 


0.00205,0.72345 


m x 


0.40883 


-0.10282 


-0.46681 


m y 


0.08899 


-0.05986 


-0.54760 


m z 


0.79605 


-0.12989 


-0.07209 


m 


0.89931 


0.17614 


0.72317 



But 4-potentials that are confined to the same local sub- 
spaces and vanish everywhere else are highly pathologi- 
cal (involving delta- or step functions). Thus, only the 
well-known special cases of local saturation NU survive 
in the continuum limit, namely, the 1-electron and the 
completely polarized, ferromagnetic case. 

2. Global noncollinear NU. The number of constraints, 
see Eq. lfT2]l. that need to be imposed on ground states 
to exhibit this type of NU, grows with the number of 
lattice sites. Thus, global noncollinear NU becomes in- 
creasingly rare for larger lattices, and is thus ruled out in 
the continuum limit, in agreement with Ref. Q. Again, 
only the 1-electron special case survives. 



These findings are reassuring for the practical appli- 
cation of SDFT to electronic structure calculations in 
atoms, molecules and solids. In the collinear case, all 
that is required in a spin-dependent Kohn-Sham calcula- 
tion is to fix two constants in the spin-up and spin-down 
channel, for example through the asymptotic behavior of 
the potentials. In the noncollinear case, a single constant 
appears to be sufficient. In situations with full spin po- 
larization, such as in half-metallic ferromagncts, the NU 
in SDFT is likely to result in the occurrence of disconti- 
nuities in the exchange-correlation potential, which will 
require further study [^,0. 

On the other hand, spin systems on small lattices are 
of great interest in the field of spintronics and quantum 
computation. From a basic point of view, for example, 
a quantum dot molecule constitutes a two-point lattice. 
The results presented in this paper will be relevant for 
the manipulation of electronic charges and spins on such 
small lattice systems, for instance in performing qubit 
operations using external magnetic fields. 
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suggesting this problem, and Giovanni Vignale and Klaus 
Capelle for valuable discussions. 
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